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ABSTRACT 


Thirteen  Douglas  Fir  beams  were  tested  in  an  investigation  of 
the  effect  of  delamination  on  the  longitudinal  strain  distribution  in 
glue  laminated  timber  beams.  An  attempt  was  made  to  determine  the 
horizontal  shear  stress  distribution  from  the  longitudinal  strain  dis¬ 
tribution.  Ten  delaminated  beams,  3"  x  6"  in  cross  section  and  6 ' -2" 
long  were  extensively  instrumented  and  strains  recorded  at  loads  up  to 
failure.  These  beams  contained  artifical  mid-depth  de laminations  of 
various  lengths  located  symmetrically  with  respect  to  midspan.  Two  types 
of  delamination  were  investigated.  In  one  type  the  wood  surfaces  at  the 
delaminated  section  were  in  contact;  whereas  in  the  other  type  the  surfaces 
were  not  in  contact . 

Beam  test  results  indicate  non-linear  longitudinal  strain  dis¬ 
tributions  exist  at  intact  sections  in  the  vicinity  of  the  delamination. 

The  extent  of  this  warping  is  dependent  on  the  ratio  of  delamination 
length  to  beam  length.  Horizontal  shear  stress  distributions,  derived 
from  the  strain  distributions  were  on  the  whole  unsatisfactory  due  to 
the  difficulty  in  accurately  defining  the  strain  distribution  across 
the  section  depth.  The  strains  measured  at  delaminated  sections  agree 
with  theoretical  values  for  low  delamination  to  beam  length  ratios.  This 
agreement  is  not  as  satisfactory  for  higher  ratios.  Strain  distributions 
indicate  no  appreciable  friction  exists  between  contacting  delamination 
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CHAPTER  I 


INTRODUCTION 

1-1  Definition  of  Delamination 

Delaminat ion  is  a  term  used  to  define  the  separation  of  wood 
surfaces  at  glued  joints.  If  the  separation  is  in  the  wood,  even  though 
very  close  to  the  glued  joint,  it  is  termed  wood  failure  or  checking. 

1-2  Causes  of  Delamination 

The  wide  spread  use  of  quality  control  in  laminating  plants 
combined  with  an  accumulation  of  experience  in  glue  laminated  timber 
fabrication  has  decreased  the  possibility  of  faulty  manufacture  of  glulam 
members.  However,  delamination  when  caused  by  faulty  manufacture  is  the 
result  of  a  defective  glued  joint.  The  most  obvious  cause  of  this  de¬ 
fective  joint  is  the  absence  of  glue  on  adjacent  laminations,  which  is 
quite  improbable  with  the  use  of  mechanical  glue  spreaders.  A  more  con¬ 
ceivable  cause  of  the  joint  defect  is  the  non-rectangularity  of  adjacent 
laminations.  The  resulting  separation  of  the  lamination  edges  produces 
a  minor  delamination,  which  may  increase  in  width  under  service  conditions. 

Shrinkage  is  a  major  cause  of  delamination.  Stress  produced 
by  deformation  may  be  built  up  as  a  result  of  this  shrinkage.  Dietz, 

Grinsf elder  and  Reissner  (1946)*  developed  a  mathematical  analysis  of 
the  stresses  produced  by  shrinkage  or  swelling  in  a  beam  composed  of 

*References  are  listed  in  the  Bibliography 
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two  laminations.  Their  analysis  indicated  that  high  stresses  can  be 
induced  by  relatively  small  changes  in  moisture  content. 

In  this  paper  two  cases  are  considered.  The  first  case  is  that 
of  a  beam  with  two  laminations,  one  of  which  is  subjected  to  a  uniform 
moisture  change.  This  moisture  change  tends  to  increase  or  decrease  the 
lamination  width  while  the  other  lamination  remains  unchanged.  As  a 
result  transverse  shear  stresses  develop  with  a  distribution  as  shown  in 
FIGURE  1.1A.  The  maximum  intensity  of  stress  at  the  edge  of  a  lamination 
is : 

&  max  =  0 . 7  e  o<  \/ExGxy 

where  e  =  o (  (M2  -  M^) 

=  coefficient  of  moisture  expansion 
M-^  =  initial  moisture  content 

M2  =  final  moisture  content 

Ex  =  modulus  of  elasticity  parallel  to  the  glue  line 
Gxy  =  shear  modulus  in  plane  of  beam  cross  section 

Using  accepted  material  properties  for  Douglas  Fir,  Dietz, 
Grinsfelder  and  Reissner  determined  the  magnitude  of  stresses  produced 
by  a  4%  change  in  moisture  content.  The  resulting  maximum  transverse 

shear  stress  developed  was  165  psi  for  edge  grain  laminations  and  240  psi 

/ 

for  flat  grain  laminations. 

The  second  case  considers  a  moisture  change  over  the  total  depth 
of  the  two  layer  beam.  In  this  case  only  normal  stresses  develop  perpen- 
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GLUE  LINE  STRESSES  RESULTING 
FROM  MOISTURE  CONTENT  VARIATIONS 

[  AFTER  DIETZ,  GRINDSFELDER  AND  REISSNER  ] 


TRANSVERSE 

shear 

STRESSES 


[  A  ]  CASE  ONE 

Cross  section  through  beam  with  one  lamination 

TENDING  TO  EXPAND  UNIFORMLY  DUE  TO  UNIFORM  INCREASE 
IN  MOISTURE  CONTENT. 


I  B  ]  CASE  TWO 

Cross  section  through  beam  when  both  laminations 

CONTRACT  DUE  TO  A  PARABOLIC  DECREASE  IN  MOISTURE  CONTENT. 


FIGURE  l.l 
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dicular  to  the  glue  line,  with  a  distribution  as  shown  in  FIGURE  1.1B. 

For  wood,  an  anisotropic  material,  the  maximum  normal  stress  at  the  edge 
of  the  lamination  is: 

CT max  =  °.45e  "<  \/Ex  ^ 

where  e  =  (M2  -  M-^) 

M2  =  moisture  content  at  the  outer  fiber 
=  moisture  content  at  the  glue  line 
Ey  =  modulus  of  elasticity  perpendicular  to  the  glue  line 

A  47o  change  in  moisture  content  results  in  a  maximum  normal  stress  of 
235  psi  for  edge  grain  laminations  and  370  psi  for  flat  grain  laminations. 

For  the  majority  of  beams  composed  of  several  laminations,  the 
usual  condition  is  intermediate  between  these  two  cases  of  moisture  change 
and  grain  orientation.  The  4%  moisture  change,  used  in  computing  normal 
and  transverse  shear  stresses,  is  quite  realistic  since  most  specifications 
commonly  allow  adjacent  laminations  to  differ  by  about  4 %  in  moisture 
content.  Canadian  Standards  Association  0122-1959  "Specification  for  Glue 
Laminated  Softwood  Structural  Timber"  allows  a  difference  of  5 %  in  moisture 
content  of  adjacent  laminations  in  a  member. 

Field  observations  indicate  that  glulam  members  tend  to  delaminate 
more  frequently  when  moisture  content  is  reduced  than  when  it  is  increased. 

1-3  Forces  in  Delaminated  Glulam  Beams 

In  an  intact  glulam  beam  the  total  compressive  or  tensile  force 
at  any  cross  section  can  be  computed  by  the  simple  flexure  formula.  This 
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force  must  be  resisted  by  the  horizontal  shearing  forces  on  the  area  of  the 
neutral  surface,  between  the  support  and  the  cross  section  in  question. 

When  a  delamination  exists  at  the  cross  section,  equilibrium  between 
flexural  forces  and  the  summation  of  horizontal  shearing  forces  must 
satisfy  equilibrium  requirements.  However,  the  computation  of  the 
flexural  force  is  made  more  difficult  due  to  the  presence  of  the  delami¬ 
nation  through  the  cross  section. 

Huggins,  Aplin  and  Palmer  (1964)  developed  a  theoretical  method 
of  analysis  to  determine  the  flexural  compressive  and  tensile  forces,  at 
any  interior  section  of  a  glulam  beam  containing  a  full  width  delamination. 
The  following  simplifying  assumptions  were  made: 

1.  A  plane  transverse  section  at  the  end  of  the  delamination 
remains  plane  after  bending. 

2.  The  radii  of  curvature  of  the  portions  of  the  beam  below  and 
above  the  delamination  are  equal. 

3.  The  moduli  of  elasticity  in  compression  and  tension  parallel 
to  the  grain  are  equal. 

4.  No  friction  exists  between  the  two  delamination  portions  of 
the  beam. 

Since  the  two  portions  of  the  delaminated  beam  have  the  same 
radius  of  curvature; 
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This  implies  that  the  bending  moments  in  the  top  and  bottom 
portions  of  the  beam  are  proportional  to  the  respective  moments  of  in¬ 
ertia  . 

For  a  beam  subjected  to  a  single  concentrated  load  as  shown  in 
FIGURE  1.2,  bending  moment  is  a  linear  function  of  vertical  shearing 
force.  Therefore,  shear  forces  in  the  top  and  bottom  portions  of  the 
beam  will  be  proportional  to  the  respective  moments  of  inertia.  From 
FIGURE  1.2  it  follows  that: 


vT 


( l-k) 3 


(l-k)3 


Mb 


vb 


If  the  total  shear  force  across  a  delaminated  section  is  known, 
Vj  and  Vg  can  be  computed. 

In  order  to  satisfy  the  condition  of  horizontal  equilibrium  the 
compressive  force  C  must  be  equal  to  the  tensile  force  T  at  a  transverse 
section  in  a  delaminated  glulam  beam.  (FIGURE  1.2).  All  quantities  in- 
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FORCES  IN  A  DELAMINATED  GLULAM  BEAM 


FIGURE  1.2 
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volved  may  be  expressed  in  terms  of  the  tensile  force  T  and  the  upper 
beam  bending  moment  M-p. 

One  equation  is  obtained  by  taking  moments  about  section  B-B 
on  the  line  of  action  of  the  force  T  as  shown  in  FIGURE  1.2. 

Rl  (a+z)  -  Mj  -  Mg  -  £4  =  0  (1.1) 

A  second  equation  is  obtained  by  considering  two  expressions 
for  the  angle  change  A0.  This  angle  change  is  the  change  in  slope  of 
the  elastic  curve  between  section  A-A  and  section  B-B  in  FIGURE  1.2. 

If  Mpx  is  the  bending  moment  in  the  top  portion  of  the  beam  at  any  point, 
a  distance  X  to  the  left  of  section  B-B,  then: 


Mtx  —  “  V>j>  X 


AG 


n 

Mrx  dx 


u 


EIt 


(1.2) 


Another  expression  for  AG  is  derived  from  a  consideration  of 
strains  in  the  top  and  bottom  portions  of  the  delaminated  beam,  under  the 
action  of  forces  C  and  T  .  Assuming  the  part  of  the  beam  to  the  left  of 
section  A-A  is  rigid,  the  change  in  length  at  mid-depth  of  the  top  portion 
of  the  beam  due  to  the  compressive  force  C  will  be: 
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Similarly,  for  the  bottom  portion  of  the  beam  the  change  in  length  due 
to  the  tensile  force  T  will  be: 


TZ 

abe 


The  change  in  length  of  the  bottom  portion  of  the  beam  relative  to  the 
top  portion  of  the  beam  is  therefore: 


CZ  +  TZ 


ATE  V 


This  relative  change  in  length  divided  by  half  the  total  beam  depth  gives 
another  expression  for  the  angle  change  AG  . 


2 

A©  =  d 


CZ  TZ 

ate  + 


abe 


(1.3) 


Combining  expressions  1.2  and  1.3  we  obtain: 


2Z 

d 


A’j-E 


aBE 


MfX  dx 

EIt 


(1.4) 


From  the  simultaneous  solution  of  Equations  1.1  and  1.4  the 
tensile  force  T  and  bending  moment  Mr  can  be  evaluated.  Subsequently, 
after  the  determination  of  bending  moment  MB  and  compressive  force  C  , 
the  distribution  of  stress  or  strain  can  be  found  at  any  delaminated 
cross  section.  With  the  compressive  force  C  known  an  average  horizontal 


shear  stress  along  the  undelaminated  length  of  the  plane,  in  which  the 
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delamination  occurs,  can  be  computed. 

1-4  Shear  Stress  Concentration  in  Delaminated  Glulam  Beams 

There  are  two  methods  of  analysis  which  may  be  applied  to 
postulate  a  stress  concentration  factor  for  the  shear  stress  near  the 
end  of  a  de lamination. 

The  first  method  is  based  on  work  by  Guyon  (1960)  related  to 
anchorage  zone  stresses  in  post -tensioned  prestressed  concrete  beams. 
Considering  local  stresses  near  the  end  surface  of  a  homogeneous,  isotropic, 
prismatic  member  produced  by  a  concentrated  load  acting  on  this  end  surface, 
Guyon  determined  influence  coefficients  for  calculating  shear  stresses 
along  the  member.  A  loading  condition  very  similar  to  the  condition 
analysed  by  Guyon  exists  in  a  vertical  plane  at  the  end  of  a  delamination 
in  a  glue  laminated  beam.  The  flexural  stress  distribution  in  the  two 
portions  of  a  delaminated  beam,  proposed  by  Huggins,  Aplin  and  Palmer, 
satisfies  this  condition  of  concentrated  loading  when  the  stress  distri¬ 
bution  is  sub-divided  into  a  series  of  concentrated  loads,  extending 
across  the  beam  depth.  (FIGURE  1.3). 

However,  Guyon 1 s  solution  for  anchorage  zone  shear  stresses  is 
only  approximate  as  indicated  by  Iyengar  (1962).  Iyengar  investigated 
the  problem  of  anchorage  zone  stresses  by  employing  a  two  dimensional 
model.  The  resulting  solution,  an  Airy  stress  function,  has  been  evalu¬ 
ated  by  Lee  (1964).  The  resulting  influence  lines  for  shear  stresses 
at  various  positions  along  the  neutral  axis  of  a  beam  are  shown  in 
FIGURE  1.4.  This  second  method  of  computing  shear  stress  concentrations 
has  the  one  disadvantage,  in  that  it  can  be  applied  only  to  the  case 
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FORCES  AND  STRESSES  AT  THE 
TERMINATION  OF  A  DELAMINATION 
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FIGURE  1.3 
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INFLUENCE  LINES  FOR  SHEAR  STRESSES 


Shear  Stress  =  CP 

2B 


FIGURE  1.4 
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where  the  delaminat ion  occurs  at  the  neutral  axis  of  the  beam.  Guyon's 
method,  on  the  other  hand,  may  be  applied  regardless  of  the  position  of 
the  delamination. 

1-5  Distribution  of  Stress  and  Strain  in  Douglas  Fir  Glulam  Beams 

The  elastic  beam  theory  predicts  with  certain  limits  the  flexu¬ 
ral  behavior  of  wood  beams.  The  theory  may  be  applied  to  wood  beams  if  the 
tensile  strength  is  reduced  by  defects,  so  that  tensile  failure  occurs  before 
the  compression  stress  reaches  the  proportional  limit  value.  At  the  critical 
section  of  the  beam  the  distribution  of  stress  and  strain  is  linear  up  to 
failure.  The  proportional  limit  stress  in  bending  for  clear  wood  is  much 
greater  than  the  value  determined  by  the  compression  test  parallel  to  the  grain. 
In  bending  it  usually  exceeds  the  crushing  strength  determined  from  the 
compression  test.  As  shown  below  a  linear  stress -strain  relation,  on 
which  the  elastic  theory  is  based,  is  not  valid  for  these  beams  when  the 
failure  load  is  approached. 

When  a  clear  wood  beam  is  subjected  to  bending  deformation, 
within  the  elastic  limit,  the  strain  distribution  and  corresponding  stress 
distribution  across  the  depth  of  the  beam  are  as  shown  in  FIGURE  1.5A. 

With  progressive  bending  deformation  the  compressive  strain  increases 
beyond  the  proportional  limit  strain.  As  a  result  the  extreme  compres¬ 
sion  fibers  begin  to  deform  in  a  non-linear  elastic  manner.  (FIGURE  1.5B). 
Further  bending  causes  the  compression  fibers  to  reach  their  elastic 
limit  and  deform  plastically  while  the  tension  fibers  continue  to  deform 
elastically.  Consequently  a  stress  redistribution  occurs  and  the  neutral 
axis  shifts  down  toward  the  tension  face  as  in  FIGURE  1.5C.  This  stress 
redistribution  continues  with  increasing  applied  moment  until  a  tension 
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GLULAM  STRESS  -  STRAIN  DISTRIBUTIONS  FOR 
PROGRESSIVE  BENDING  DEFORMATIONS 


t  A  ]  BENDING  WITHIN  THE  PROPORTIONAL  LIMIT 
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[  B  1  BENDING  AT  THE  PROPORTIONAL  LIMIT 


[  A  ]  BENDING  IN  THE  PLASTIC  REGION 


FIGURE  1. 5 
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failure  occurs. 

Ramos  (1961)  investigated  strain  distribution  in  Douglas  Fir 
beams  within  the  plastic  range.  FIGURE  1.6  illustrates  the  strain  dis¬ 
tributions  obtained  preceding  failure  for  three  test  beams  2"  x  6"  in 
cross  section  with  a  span  of  12' -0".  From  this  work  it  was  concluded 
that  the  assumption  of  a  linear  distribution  of  longitudinal  strain  is 
valid  for  practical  purposes.  The  strain  distributions  in  FIGURE  1.6 
clearly  illustrate  the  gradual  shift  of  the  neutral  axis  towards  the 
tension  face  of  the  beam. 

1-6  Previous  Experimental  Investigations 

One  of  the  major  problems  in  the  area  of  timber  engineering  lies 
in  the  evaluation  of  the  strength  characteristics  of  beams.  This  is  very 
evident  when  studies  on  the  horizontal  shear  strength  of  glulam  beams 
are  examined  for  conclusive  results.  Furthermore,  there  are  a  limited 
number  of  published  investigations  concerning  horizontal  shear  stresses 
in  delaminated  glulam  beams. 

Huggins,  Aplin  and  Palmer  (1964)  tested  175  beams,  3  1/4"  x 
6"  x  9'-0",  under  static  and  dynamic  concentrated  loads  at  the  centerline 
of  a  8'-0"  span.  One  hundred  and  fifteen  of  these  beams  contained  full 
width  delaminations  at  various  locations  in  the  beam.  It  was  observed 
that  an  average  shear  stress,  calculated  on  the  undelaminated  surface, 
did  not  explain  the  relative  incidence  of  horizontal  shear  failures  for 
the  different  delamination  patterns.  Shear  stress  concentrations  calcu¬ 
lated  by  Guyon's  method  were  considerably  lower  than  those  obtained  by 
Lee's  method.  Using  this  latter  method  very  good  correlation  was  ob- 
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LONGITUDINAL  STRAIN  DISTRIBUTIONS 

[  AFTER  RAMOs] 
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PLANE  2 


NOTE 

PLANES  I  AND  2  REFER  TO  TRANSVERSE  SECTIONS 
IN  THE  CONSTANT  MOMENT  REGION  OF  THE  BEAMS 


FIGURE  1.6 
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tained  between  the  shear  strength  of  glulam  beams,  indicated  by  shear 
block  tests,  and  the  occurrence  of  horizontal  shear  failures.  The  re¬ 
ducat  ion  of  strength  under  repeated  loading  was  confirmed.  However, 
the  extent  of  this  occurrence  was  impossible  to  determine,  due  to  the 
manner  of  the  test  and  the  limited  number  of  beams  tested. 

A  previous  investigation  at  the  University  of  Alberta,  con¬ 
cerning  the  effect  of  delamination  on  the  shear  strength  of  glulam  beams, 
was  conducted  by  Strong  (1964).  In  this  comprehensive  program,  239 
Douglas  Fir  beams  and  110  shear  block  specimens  were  tested.  Length  of 
delamination  and  position  of  delamination  were  the  major  variables  in 
the  3"  x  6"  x  6 ' -2"  long  beams.  Artificial  delaminations  were  introduced 
by  inserting  two  layers  of  wax  paper  at  the  glue  line.  Test  results 
indicated  that  an  increase  in  length  of  delamination  had  the  effect  of 
decreasing  the  load  carrying  capacity  of  the  beams.  A  de lamination 
located  at  mid-depth  was  shown  to  be  more  critical  than  a  delamination 
at  1/4  depth  or  3/4  depth.  End  delamination  proved  to  be  more  critical 
than  centerline  delamination.  As  in  previous  investigations  no  corre¬ 
lation  was  found  to  exist  between  shear  block  and  beam  tests. 

In  the  previous  experimental  investigations  described  horizontal 
shear  stresses  in  delaminated  glulam  beams  have  been  in  terms  of  the 
average  stress  sustained  prior  to  failure.  However,  if  shear  stress 
concentrations  do  exist  at  the  termination  of  a  delamination,  they  must 
have  an  effect  on  longitudinal  strain  distributions  in  this  vicinity. 
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1-7  Scope  of  Present  Investigation 

The  investigation  described  in  this  thesis  was  conducted 
primarily  for  the  purpose  of  studying  the  effect  of  delamination  on  the 
longitudinal  strain  distribution  in  glue  laminated  Douglas  Fir  beams. 
Secondly,  the  feasibility  of  calculating  the  actual  horizontal  shear 
stress  distribution  from  the  strain  distribution  was  examined. 

Thirteen  rectangular  simply  supported  beams  were  tested  to 
failure  under  the  application  of  a  single  concentrated  load  at  midspan. 
Ten  of  the  beams  contained  artificial  delaminations  in  which  the  type 
of  delamination  and  length  of  delamination  were  major  variables.  In 
these  beams  the  location  of  the  delamination  was  at  mid-depth  of  the 
section  and  symmetrical  with  respect  to  the  centerline  of  the  span. 
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CHAPTER  II 


HORIZONTAL  SHEAR  STRESS  DISTRIBUTION 

2-1  Introduction 

In  the  development  of  the  expression  for  flexural  stresses  in 
beams  of  homogeneous  material  it  is  assumed  that  transverse  sections 
which  are  plane  before  bending  remain  plane  during  bending.  For  a 
material  which  obeys  Hooke's  Law  this  results  in  the  familiar  equation: 

f  =  MX 
I 


Since  horizontal  shearing  stresses  are  related  to  the  distribution  of 
flexural  stresses,  it  follows  that  the  same  assumption  of  plane  sections 
remaining  plane  is  made  in  deriving  expressions  For  shear  in  beams  of 
homogeneous  material.  In  other  words  the  equation  for  shearing  stress 
in  a  rectangular  beam 


v 


is  based  on  the  assumption  that  sections  do  not  warp  due  to  bending. 
However,  if  warping  of  the  section  occurs,  the  horizontal  shear  stresses 

VQ 

cannot  be  expressed  as  rrr  . 
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2-2  Warped  Transverse  Sections 

Consider  a  simple  span  rectangular  beam  with  a  single  con¬ 
centrated  load  at  midspan,  as  shown  in  FIGURE  2.1.  If  we  assume  that 
due  to  bending  transverse  sections  of  the  beam  are  warped,  we  may  ex¬ 
press  the  longitudinal  strains  by  the  equation: 

l 

e  =  ay  n 

where  y  =  the  distance  from  the  neutral  axis  to  a  longitudinal 

fiber 

e  =  the  longitudinal  strain  at  this  fiber 
a  =  a  constant 

n  =  an  index  which  determines  the  degree  of  the  equation 
defining  the  strain  distribution 

Assuming  stress  is  proportional  to  strain,  the  change  in  flexural  stress 
on  any  fiber  between  section  1  and  2  (FIGURE  2. 1C)  is: 

1 

Af  =  -  a2)y  n 

1 


Considering  the  horizontal  equilibrium  of  the  element  ABCD,  we  can  derive 
the  expression  for  horizontal  shear  stresses  v  . 

v  bAX  = 
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(2.1) 


If  ft  =  1  in  the  above  expression  the  horizontal  shear  stress 
becomes  equal  to  the  familiar  ~  .  This  distribution,  shown  in 

FIGURE  2.2  has  an  ordinate  at  the  neutral  axis  equal  to  1.5  times  the 
average  shear  stress  ^  . 

To  illustrate  the  effect  of  warping  let  us  now  take  n  =  2  . 
To  satisfy  the  equilibrium  of  internal  and  external  moments: 


where 


e 

L 


=  outer  fiber  strain  for  a  parabolic  strain  distribution 


=  outer  fiber  strain  for  a  linear  strain  distribution 
e  =  ay 


Using  this  expression  for  ep  together  with  equation  2.1  the  horizontal 
shear  stress  distribution  across  the  section  can  be  computed  and  is  shown 
in  FIGURE  2.2. 

FIGURE  2.2  indicates  that  warping  of  cross  sections  increases 
the  magnitude  of  the  maximum  horizontal  shear  stress.  Therefore,  if 
longitudinal  strains  can  be  measured  and  the  strain  distribution 
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HORIZONTAL  SHEAR  STRESS  DISTRIBUTIONS 
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established,  it  should  be  possible  to  compute  the  actual  shear  stress 
distributions . 


CHAPTER  III 


MATERIAL  PROPERTIES 

3-1  Moisture  Content  and  Specific  Gravity 

Prior  to  the  fabrication  of  the  beam  specimens  and  tension  and 
compression  test  specimens  all  laminating  stock  was  checked  for  exces¬ 
sive  variations  in  moisture  content  by  means  of  an  electrical  resistance 
moisture  indicator.  The  "Moore  Moisture  Meter  Model  MN-1"  used  for  this 
purpose  is  shown  in  PLATE  3. LA.  The  moisture  content  of  the  laminating 
stock  used  in  this  program  was  controlled  between  10%  and  12%  with  an 
average  of  11.5%. 

After  fabrication  the  beam  specimens  were  wrapped  in  polyethyl¬ 
ene  to  maintain  the  beam  moisture  content.  On  completion  of  a  beam  test 
one  moisture  content  and  specific  gravity  sample  ,  measuring  about 
3/4"  x  3"  x  3"  ,  was  cut  from  each  lamination  through  the  depth  of  the 
beam.  The  samples  were  taken  from  the  region  of  the  beam  in  which  failure 
was  initiated.  The  moisture  content  and  specific  gravity  of  each  sample 
were  determined  by  methods  outlined  in  APPENDIX  A. 

3-2  Modulus  of  Elasticity  in  Tension 

A  total  of  five  tension  test  specimens,  selected  from  flat 
grained  fabricating  stock,  were  tested  to  determine  the  modulus  of 
elasticity  parallel  to  the  grain  in  tension.  FIGURE  3.1  indicates 
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28. 


the  nominal  dimensions  for  these  tension  test  specimens. 

All  specimens  were  tested  in  a  "Baldwin-Tate-Emery  Universal 
Testing  Machine"  of  200,000  pound  capacity.  Strains  were  measured  by 
means  of  a  "Model  TS-M  Dual  Extensometer"  (PLATE  3. IB),  with  an  8-inch 
gauge  length,  in  conjunction  with  an  automatic  electrical  load -deforma¬ 
tion  recorder.  At  about  80%  of  the  ultimate  load  slippage  of  the  speci¬ 
men  in  the  testing  machine  grips  occurred.  As  valid  load  deformation 
plots  were  impossible  to  obtain  at  this  stage,  the  extensometer  was  re¬ 
moved  and  loading  continued  to  failure.  All  specimens  exhibited  a 
splintering  tension  type  of  failure.  After  failure,  moisture  content 
and  specific  gravity  determinations  were  performed. 

3-3  Modulus  of  Elasticity  in  Compression 

Five  compression  test  specimens,  measuring  2"  x  2"  x  l'-2", 
were  tested  to  determine  the  modulus  of  elasticity  in  compression 
parallel  to  the  grain.  These  specimens  were  cut  from  the  end  regions 
of  XY-4  beam  specimens  previously  tested  to  failure.  Longitudinal  strains 
over  an  8-inch  gauge  length  were  obtained  by  means  of  an  "H.F.  Moore 
Mechanical  Extensometer".  After  crushing  failure  occured,  moisture  con¬ 
tent  and  specific  gravity  values  were  obtained  for  all  compression  test 


specimens . 
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CHAPTER  IV 


BEAM  SPECIMENS  AND  FABRICATION 

4-1  Beam  Specimen  Types 

Thirteen  beams,  3"  x  6"  in  cross  section  and  6' -2"  long,  were 
tested  in  this  program.  These  beams  were  sub-divided  into  three  groups 
as  described  in  TABLE  4-1 .  The  first  group,  consisting  of  three  control 
beams,  contained  no  de laminations .  The  ten  delaminated  beams  were 
divided  into  two  groups  each  comprising  five  different  types.  In  one 
group  of  delaminated  beams,  the  wood  surfaces  at  the  delaminated  section 
were  in  contact;  whereas  in  the  other  group  the  surfaces  were  not  in 
contact.  Details  are  shown  in  FIGURE  4.1. 

4-2  Beam  Fabrication  Procedures 

Prior  to  beam  fabrication  the  laminating  stock  was  examined  and 
sorted  with  respect  to  moisture  content,  approximate  specific  gravity, 
grain  orientation  and  occurrence  or  knots  and  checks.  The  resulting  beams 
consisted  of  eight  3/4"  clear,  check  free  laminations  of  Grade  B  Douglas 
Fir  laminating  stock.  These  laminations  were  generally  flat  grained. 

The  moisture  content  was  controlled  to  a  tolerance  of  +  1%  and  the  specific 
gravity  varied  from  a  maximum  at  beam  mid-depth  to  a  minimum  at  the  outer 
fiber  s . 

In  order  to  facilitate  convenient  manufacture  of  the  delaminated 
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beams,  which  contained  an  actual  clearance  between  laminations,  a  special 
fabrication  procedure  was  followed.  The  top  and  bottom  halves  of  these 
beams,  each  consisting  of  four  laminations,  were  individually  glued  and 
clamped  in  the  jigs  as  separate  units.  (PLATE  4. 1A)  .  After  allowing 
twenty  four  hours  for  curing  of  the  casein  glue,  the  beam  halves  were 
removed  from  the  jigs  and  the  delamination  prepared.  On  each  beam  half, 
material  to  a  depth  of  1/16  inch  was  removed  from  the  lamination  adjacent 
to  the  future  mid-depth  glue  line  of  the  beam  specimen.  This  operation 
was  performed  along  the  full  length  of  the  required  delamination  by  a 
power  router  having  a  cutting  head  diameter  of  6  inches.  This  preparation 
resulted  in  a  very  clean  full  width  gouge  of  uniform  depth  having  an  end 
radius  of  3  inches.  The  beam  halves  were  glued  on  the  unrouted  surfaces 
of  the  mid-depth  laminations,  assembled  as  full  beams  in  the  jigs  and 
cured  for  twenty  four  hours. 

The  second  group  of  delaminated  beams,  containing  the  wood  con¬ 
tact  type  delamination,  were  fabricated  following  the  same  basic  pro¬ 
cedure  used  for  the  clearance  type  delamination.  However,  instead  of 
routing  out  material  to  form  the  artificial  delamination  the  following 
method  was  substituted.  Initially  3  inch  square  double  layers  of  wax 
paper  were  stapled  at  the  ends  of  the  required  delamination.  After 
glue  was  spread  from  each  end  of  the  beam  halves,  up  to  and  not  including 
the  wax  paper,  the  halves  were  assembled.  PLATE  4.  IB  shows  the  glued 
beam  halves  prior  to  full  beam  assembly  and  clamping  in  the  jigs. 

To  ensure  uniformity  of  fabrication,  the  two  stage  gluing  process 
was  also  employed  in  fabricating  the  control  beams.  However,  at  the 
second  gluing  stage  the  entire  surfaces  of  the  beam  halves  were  glued 
resulting  in  a  continuous  glue  joint  at  mid-depth  of  the  control  beams. 
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[A  ] 


BEAMS  IN  CLAMPING  JIGS 


[  B  ]  HALF  BEAMS  AFTER  GLUING 


PLATE  4.1 
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CHAPTER  V 


BEAM  TESTS 


5-1  Loading  System 

The  beam  testing  program  was  conducted  using  a  standard  testing 
frame  and  the  "Amsler  Hydraulic  Jack  Plant"  in  the  Graduate  Structures 
Laboratory.  FIGURE  5.1  illustrates  the  test  equipment  arrangement. 

Load  was  applied  at  the  midspan  through  a  3"  x  4"  x  3/8"  dis¬ 
tribution  plate.  When  testing  the  delaminated  beams,  containing  the 
clearance  type  delamination,  3"  x  4"  lubricated  sheet  metal  shims  were 
snugly  fitted  in  the  delamination  at  midspan.  This  prevented  local  bend¬ 
ing  in  the  top  half  of  the  beam  and  assured  equal  load  distribution  to 
both  halves  of  the  beam. 

Strong,  in  the  previous  investigation  at  the  University  of  Alberta, 
encountered  difficulties  in  testing  a  number  of  beams  to  failure  because 
of  lateral  instability.  In  the  present  project  reaction  blocks,  as 
shown  in  PLATE  5.1A,  were  provided  to  prevent  lateral  buckling.  These 
blocks,  one  hinged  and  one  mounted  on  a  roller,  satisfied  fully  the  con¬ 
ditions  of  simple  supports  while  still  providing  lateral  support  to  the 
beam.  The  reactions  were  transferred  through  these  reaction  blocks  to 
the  beam  through  2"  x  3"  x  3/8"  bearing  plates. 
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BEAM  TESTING  ARRANGEMENT 


FIGURE  5.1 
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[  A  ]  TYPICAL  BEAM  TESTING  ARRANGEMENT 


[  B  ]  CONVERTED  DEMEC  GAUGE 


PLATE  5.1 


1  t&AT  MA38  JAD49YT  [  A  ] 
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5-2  Beam  Instrumentation 

i .  Deflection  Measurement 

Vertical  deflections  were  measured  at  one  foot  intervals 
along  the  beam  by  dial  gauges  having  a  one-inch  travel  and  reading  to 
0.001  inch. 

ii .  Demec  Strain  Measurements 

In  order  to  measure  strains  over  a  small  gauge  length  it  was 
necessary  to  modify  an  8-inch  Demec  strain  gauge  to  obtain  readings  on 
a  2 -inch  gauge  length.  This  was  achieved  by  clamping  a  removable 
conversion  bar  securely  to  the  Demec  gauge  as  shown  in  PLATE  5. IB. 

This  bar  contained  a  machined  hole  at  one  end  which  fitted  over  the  fixed 
Demec  gauge  point,  assuring  no  relative  movement  between  the  gauge  and 
conversion  bar.  At  6 -inches  from  the  fixed  Demec  gauge  point,  another 
fixed  gauge  point  was  located  on  the  conversion  bar.  This  modification 
made  it  possible  to  record  deformations  over  a  2-inch  gauge  length  with 
an  accuracy  of  +  0.000004  in. /in.  A  calibration  factor  for  the  2-inch 
Demec  strain  gauge  was  obtained  from  a  tension  test  of  a  1/2  inch  diameter 
steel  rod.  This  rod  was  instrumented  with  two  electrical  resistance  SR-4 
Type  A -3  strain  gauges  and  two  sets  of  Demec  gauge  points,  providing  the 
necessary  data  for  the  computation  of  a  calibration  factor. 

Demec  gauge  points  were  bonded  with  sealing  wax  to  the  dela¬ 
minated  beams  and  control  beam  X-003  in  locations  as  shown  in  FIGURE  5.2 


and  5.3  respectively. 
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iii .  Electrical  Resistance  Strain  Measurement 

In  order  to  provide  a  check  on  the  2 -inch  Demec  strain  gauge 
readings,  electrical  resistance  SR-4  Type  A-3  strain  gauges  were  employed 
on  beams  X2-41  and  X2-42. 

In  the  region  of  electrical  strain  gauge  instrumentation 
(FIGURE  5.2)  the  beam  surface  was  sanded  and  shellac  was  applied  in  pre¬ 
paration  for  mounting  the  SR-4  gauges.  The  shellac  acted  as  a  sealer  in 
protecting  the  sensitive  SR-4  gauges  from  moisture  in  the  wood.  The  SR-4 
Type  A-3  gauges  were  then  attached  to  the  beam  using  the  bonding  and  curing 
compounds  supplied  in  the  "Budd  GA-1  Cement  Kit"  according  to  the  manu¬ 
facturer's  recommended  procedure.  All  gauges  were  checked  with  an 
ohmeter  before  and  after  wiring  and  all  defective  gauges  were  replaced. 

A  problem  of  drift  appeared  on  attempting  to  read  strains  using 
a  "Baldwin -Lima -Hamilton  20  Point  Switching  Unit"  in  conjunction  with  a 
"Baldwin -Lima -Hamilton  Type  N,  SR-4  Strain  Indicator".  An  initially 
balanced  reading  on  the  SR-4  indicator  began  to  drift  immediately  after 
the  particular  SR-4  gauge  was  switched  into  the  bridge  circuit.  FIGURE 
5.4  illustrates  the  manner  in  which  the  strain  readings  drifted  with 
time.  This  problem  of  gauge  drift  can  be  solved  two  ways.  The  first 
is  to  power  the  gauge  continuously  avoiding  the  drift  with  each  reading. 

The  second  method  is  to  take  strain  readings  at  a  constant  time  interval 
after  the  gauge  has  been  powered.  This  latter  method  of  obtaining  a 
stabilized  reading  was  employed  in  this  investigation.  The  strain  in¬ 
dicator  was  balanced  and  strains  were  recorded  after  a  two  minute  warm 


up  for  each  strain  gauge. 
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5-3  Testing  Procedure 

Prior  to  the  application  of  load  the  initial  Demec  strain  gauge 
readings  were  recorded  from  the  end  of  the  beam  to  the  vertical  center- 
line  of  the  beam.  Where  applicable,  the  SR-4  gauges  were  read  following 
the  reading  of  the  Demec  gauge  at  the  corresponding  section.  An  initial 
load  of  150  pounds  was  applied  before  the  dial  gauges  were  read.  Loading 
was  then  applied  at  a  rate  of  250  pounds  per  minute  until  the  load  incre¬ 
ment  of  1000  or  1500  pounds,  depending  on  the  beam  type,  was  obtained. 

On  the  delaminated  beams  a  small  amount  of  sealing  wax  was 
placed  at  the  termination  of  the  artificial  delamination  to  indicate 
readily  any  separation  of  the  glue  joint  and  further  propagation  of 
the  delamination.  At  each  load  increment  the  ends  of  the  delamination 
were  examined  for  cracking  before  strain  and  dial  gauges  were  read. 

If  at  any  stage  cracking  occurred  in  the  sealing  wax,  recording  of  strain 
was  terminated  and  the  loading  was  increased  to  failure. 

In  the  later  stages  of  loading  the  usual  non-linear  relation 
between  load  and  deflection  developed.  This  necessitated  a  waiting  period 
of  5  to  10  minutes  before  recording  of  strains  could  commence.  Occasional 
checks  on  strains  in  regions  of  maximum  moment,  before  and  after  this 
waiting  period,  showed  an  almost  imperceptible  change  in  the  outer  fiber 
strains.  The  duration  of  a  beam  test  varied  from  4  to  6  hours. 
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CHAPTER  VI 


PRESENTATION  OF  TEST  RESULTS 

6-1  Tension  and  Compression  Test  Results 

The  test  results  of  the  laminating  stock  tension  and  compres¬ 
sion  test  specimens  are  presented  in  FIGURES  6.1,  6.2  and  TABLE  VI-1. 

6-2  Beam  Test  Results 

Beam  test  results  are  presented  in  tabular  form  as  follows: 

1.  Control  beam  test  data  (TABLES  VI-2  and  VI-3) 

2.  Delaminated  beam  test  data  (TABLES  VI-4  and  VI-5) 

3.  Computed  shear  stresses  (TABLE  VI- 6) 

Further  test  results  are  presented  in  graphical  form  as  follows: 

1.  Control  beam  strain  distributions  (FIGURE  6.3) 

2.  Delaminated  beam  strain  distributions  (FIGURES  6.4  to  6.20) 

3.  Load-deflection  curves  (FIGURES  6.21  to  6.24) 

4.  Horizontal  shear  stress  distributions  (FIGURES  6.25  and  6.26). 
PLATE  6 . 1A  shows  the  typical  horizontal  shear  failure  which 

occurred  in  all  the  delaminated  beams  and  in  one  control  beam.  PLATE  6 . IB 
shows  the  typical  occurrence  of  the  horizontal  shear  failure  beginning  at 
the  glue  line  near  the  end  of  the  delamination.  Proceedings  towards  the 
reaction  the  failure  surface  propagates  into  the  lamination  rather  than 
continuing  along  the  glue  line. 
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A  computer  program  (APPENDIX  B)  was  written  in  order  to  obtain 
shear  stresses  from  horizontal  equilibrium  considerations,  using  the 
strain  distributions  obtained  in  the  tests. 
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STRESS  KSI 


STRAIN  [  in. /in.  ] 

FIGURE  6.1 


STRESS  KSI 
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DISTRIBUTIONS  FOR  A  3'-0  LONG  DELAMINATION 


SECTION  (p  (*) 


Beam  XY-43  @  P  =  5000.  lbs. 
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EXPERIMENTAL  AND  THEORETICAL  STRAIN 
DISTRIBUTIONS  FOR  A  4'-0  LONG  DELAMINATION 


Beam  XY-44  @  P  =  4000.  lbs. 
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FIGURE  6.18 
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COMPARATIVE  DEMEC  AND  SR-4 
STRAIN  DISTRIBUTIONS 


SECTION 


Beam  XZ-41  @  P-7000,  lbs. 
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Beam  XZ-42@  P-7500.  lbs. 


FIGURE  6.20 
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[  A  ]  TYPICIAL  HORIZONTAL  SHEAR  FAILURE 


[  B  ]  HORIZONTAL  SHEAR  FAILURE  SURFACE 


PLATE  6.1 


— 


CHAPTER  VII 


DISCUSSION  OF  TEST  RESULTS 

7-1  Material  Properties 

Examination  of  the  tension  and  compression  test  results  (FIGURES 
6.1,  6.2  and  TABLE  VI-l)  indicate  a  slight  difference  exists  between  the 
average  compressive  and  tensile  moduli  of  elasticity  parallel  to  grain. 
Tests  performed  at  the  United  States  Forest  Products  Laboratory  have  also 
showed  a  difference  in  the  tensile  and  compressive  moduli  of  elasticity. 
However,  this  difference  was  insignificant  for  Douglas  Fir. 

6 

The  value  of  the  average  modulus  of  elasticity  (1.72  x  10  psi) 

for  the  control  beams,  computed  from  the  initial  tangents  to  the  load- 

deflection  curves,  is  lower  than  the  moduli  of  elasticity  determined 

from  the  tension  and  compression  tests.  Examination  of  FIGURES  6.14  to 

6,18  illustrating  theoretical  and  experimental  strain  distributions  at 

£ 

delaminated  sections  imply  a  modulus  of  elasticity  between  1.7  x  10  psi 

6 

and  2.1  x  10  psi  is  appropriate. 

The  moisture  content  of  the  beam  laminations  determined  upon 
completion  of  the  beam  tests  agree  exceptionally  well  with  moisture  content 
values  obtained  prior  to  fabrication.  The  specific  gravity  values  obtained 
after  tests  exhibit  a  definite  tendency  towards  a  maximum  specific  gravity 
at  the  mid-depth  laminations,  decreasing  to  a  minimum  at  the  outer  lami¬ 
nations.  However,  a  limited  number  of  beams  did  contain  dense  material 
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in  the  outer  laminations. 

7-2  Control  Beams 


In  order  to  establish  whether  or  not  the  two  stage  method  of 
beam  fabrication  had  a  detrimental  effect  on  beam  shear  strength,  the 
control  beam  test  results  of  this  investigation  and  Strong's  investi¬ 
gation  are  compared.  In  both  these  investigations  beams  of  identical 
dimensions  were  manufactured  from  the  same  fabricating  stock.  However, 
the  beams  tested  by  Strong  were  manufactured  by  normal  fabricating  pro¬ 
cedures.  The  thirteen  control  beams  in  Strong's  investigation  sustained 
an  average  horizontal  shear  stress  of  391  psi  prior  to  failure  by  four 
typical  failure  modes.  The  three  beams  failing  in  horizontal  shear 
failed  at  an  average  horizontal  shear  stress  of  377  psi.  Examination 
of  TABLE  VI-2  indicates  all  the  control  beams  in  this  present  investi¬ 
gation  developed  horizontal  shear  stresses  comparable  to  Strong's  beam 
results.  Therefore,  the  two  stage  method  of  beam  manufacture  did  not 
have  an  apparent  effect  on  horizontal  shear  strength. 

Qualitatively  the  strain  distributions  of  control  beam  X-003 
(FIGURE  6.3)  appear  to  satisfy  the  normal  assumption  of  a  linear  dis¬ 
tribution  of  strains.  This  is  more  evident  when  a  comparison  is  made 
with  the  distributions  obtained  by  Ramos  (FIGURE  1.6).  The  strain  dis¬ 
tributions  obtained  from  control  beam  X-003  do  not  contain  the  abrupt 
variations  in  strain  through  the  depth  of  the  section,  which  are  present 
at  a  number  of  the  sections  in  the  beams  tested  by  Ramos. 

At  all  sections  in  control  beam  X-003  the  outer  tensile  strains 
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are  about  ten  percent  larger  than  the  corresponding  compressive  strains. 
Therefore,  in  order  to  satisfy  equilibrium,  this  difference  would  imply 
that  the  compressive  elastic  modulus  is  slightly  larger  than  the  tensile 
elastic  modulus. 

At  sections  1  to  4  of  control  beam  X-003,  the  neutral  axis 
shifted  upwards  from  the  centroidal  axis  of  the  beam.  Examining  TABLE 
VI- 3  with  respect  to  the  specific  gravity  of  the  laminations  in  this 
control  beam,  indicates  laminations  above  the  centroidal  axis  have  an 
average  specific  gravity  (0.579)  higher  than  the  laminations  below 
(0.522).  Doyle,  McBurney  and  Drow  (1946)  found  that  an  increase  in 
specific  gravity  resulted  in  an  increase  in  the  modulus  of  elasticity. 

A  transformed  section  concept  based  on  a  varying  modulus  of  elasticity 
would  explain  the  upward  shift  of  the  neutral  axis. 

A  computer  program  (APPENDIX  B)  was  used  in  an  attempt  to 
determine  the  horizontal  shear  stress  from  longitudinal  strain  distri¬ 
butions  in  control  beam  X-003.  This  analysis,  considering  equilibrium 
of  the  top  and  bottom  portions  of  the  beam,  usually  produces  two  incom¬ 
patible  values  of  shear  stress  at  the  neutral  axis.  Furthermore,  for  all 
loads  the  shear  stresses  at  the  neutral  axis  of  the  instrumented  sections 
generally  deviate  widely  from  theoretical  values,  as  shown  in  TABLE  VI-6B. 

7-3  Delaminated  Beams 

i ,  Strains  and  Shear  Stresses  at  Intact  Sections 

In  a  homogeneous,  isotropic,  rectangular  beam  a  linear  dis¬ 
tribution  of  longitudinal  strain  is  valid  at  the  transverse  sections 
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along  the  beam  length.  In  order  to  satisfy  equilibrium  at  any  transverse 
section,  the  shear  stress  distribution  is  a  second  degree  parabola  having 
a  maximum  ordinate  at  the  neutral  axis.  This  maximum  shear  stress  is  1.5 
times  the  average  transverse  shear  stress. 

The  presence  of  a  discontinuity  at  the  mid-depth  of  a  beam, 
such  as  a  delamination,  results  in  a  shear  stress  concentration  in  the 
vicinity  of  this  discontinuity.  To  simultaneously  satisfy  this  condition 
of  a  shear  stress  concentration  and  equilibrium  of  the  cross  section, 
a  non-linear  distribution  of  flexural  stresses  must  exist  across  the 
section.  If  the  outer  fiber  flexural  stress  is  below  the  proportional 
limit  stress,  the  strain  distribution  will  be  in  direct  proportion  to  the 
flexural  stress  distribution.  Therefore,  sections  which  are  plane 
before  bending  do  not  remain  plane,  but  warp  as  a  result  of  bending. 

Under  this  condition  the  maximum  shear  stress  at  the  neutral  axis  will 
be  in  excess  of  1.5  times  the  average  shear  stress,  as  shown  in  FIGURE  2.2. 

FIGURES  6.4  to  6.13  illustrate  the  longitudinal  strain  distri¬ 
butions  measured  in  the  instrumented  delaminated  beams.  From  these  dis¬ 
tributions  it  is  very  evident,  by  the  extent  of  the  non-linearity  in 
strains,  or  warping  of  sections,  that  a  definite  shear  stress  concentration 
must  exist  in  the  vicinity  of  the  delamination. 

FIGURES  6.4  and  6.9  indicate  parabolic  strain  distributions 
exist  at  section  4  in  beams  XZ-41  and  XY-41,  which  contain  l'-O"  long 
delaminations.  As  the  distance  from  the  end  of  the  delamination  increases 
the  strains  approach  a  linear  distribution.  Beams  XZ-42  and  XY-42 
(FIGURES  6.5  and  6.10),  containing  2 ' -0"  long  delaminat ions ,  exhibit 
marked  parabolic  strain  distributions  at  sections  4  to  6.  Similarily, 
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the  extent  and  magnitude  of  the  parabolic  strain  distributions  are  shown 
for  beams  XZ-425  and  XY-425,  containing  2 '-6"  long  delaminations ,  in 
FIGURES  6.6  and  6.11.  In  both  of  these  beams  section  4  at  the  end 
of  the  delamination  is  warped  to  such  an  extent  that  magnitude  of  the 
strains  in  the  outer  two  laminations  approach  a  constant  value  across 
these  laminations.  As  a  consequence,  a  larger  amount  of  the  total  trans¬ 
verse  shear  at  the  section  will  be  carried  by  the  central  portion  of  the 
beam.  In  beams  XZ-43  and  XY-43,  with  3'-0"  long  delaminations  (FIGURES 
6.7  and  6.12),  parabolic  strain  distributions  exist  at  sections  4  to  8. 
FIGURE  6.19B,  typical  for  all  delaminated  beams,  indicates  the  fairly 
uniform  progression  of  strains  at  an  intact  section,  through  the  full 
range  loads  to  failure.  The  strain  distributions  in  beams  XZ-44  and  XY-44, 
containing  4'-0"  long  delaminations  (FIGURES  6.8  and  6.13)  are  peculiar. 

At  sections  4  and  5  the  strain  increases  from  the  neutral  axis  to  a 
maximum  at  the  quarter  depth  before  decreasing  towards  the  outer  lami¬ 
nations.  In  the  delaminated  portions  of  these  beams  the  compression  and 
tension  flexural  forces  are  largely  concentrated  near  the  delamination 
plane,  as  shown  by  the  strains  in  FIGURE  6.18.  This  immediate  proximity 
of  the  larger  portion  of  these  longitudinal  forces  to  the  neutral  surface 
will  result  in  substantial  horizontal  shear  stress  concentrations  at 
intact  sections,  in  the  vicinity  of  the  delamination.  This  would  then 
explain  the  extensive  warping  observed  in  these  beams. 

It  is  therefore  evident  that  the  longer  the  delamination  or 
larger  the  ratio  of  delamination  length  to  beam  length  (D/L  ratio)  the 
more  the  strain  distribution,  in  the  vicinity  of  the  delamination, 
deviates  from  a  linear  distribution.  In  beams  with  low  D/L  ratios  the 
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longitudinal  strain  is  a  maximum  at  the  outer  fibers,  decreasing  almost 
linearly  to  zero  at  the  neutral  axis.  For  larger  D/L  ratios  the  strains 
at  intact  sections  in  the  vicinity  of  the  delamination  become  parabolic. 

For  these  large  D/L  ratios  the  strains  in  delaminated  sections  at  the 
end  of  the  delamination  are  a  minimum  at  the  outer  fibers,  increasing 
to  a  maximum  at  the  delamination.  Guyon's  and  Lee's  influence  lines 
(Figure  1.4)  indicate  loads  at  lower  eccentricities  produce  larger  shear 
stresses  on  the  neutral  surface.  Therefore,  beams  with  high  D/L  ratios 
will  develop  large  shear  stress  concentrations.  This  explains  the  pro¬ 
nounced  warping  observed  at  these  cross-sections.  Similarly,  beams  with 
low  D/L  ratios  having  the  larger  portion  of  the  longitudinal  fiber  forces 
at  large  eccentricities,  will  develop  less  severe  shear  stress  concentrations 
resulting  in  moderate  warping. 

Using  the  computer  program  in  APPENDIX  B,  an  attempt  was  made 
to  compute  the  horizontal  shear  stress  distributions  between  transverse 
beam  sections,  for  all  loads  to  failure.  In  the  analysis  of  every 
section  a  second  degree  parabola  is  fitted  to  each  set  of  five  recorded 
longitudinal  strains  above  and  below  the  neutral  axis.  The  flexural 
stresses  at  each  section  are  determined  from  the  defined  strains  together 
with  a  modulus  of  elasticity  having  a  value  equal  to  the  average  of  the 
compressive  and  tensile  moduli.  The  shear  stress  distribution,  between 
any  two  successive  sections,  is  then  obtained  by  considering  the  hori¬ 
zontal  equilibrium  of  thirty  elements  through  the  depth  of  the  top  and 
bottom  halves  of  the  beam.  For  the  vast  majority  of  the  sections  in 
all  the  beams  this  analysis  usually  produces  two  different  values  of 
shear  stress  at  the  neutral  axis.  (TABLE  VI— 6) .  FIGURE  6.25  indicates 
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the  horizontal  shear  stress  distributions  at  sections  where  compatabi lit y 
of  shear  stress  at  the  neutral  axis  is  approached.  These  distributions 
when  compared  to  the  normal  parabolic  distribution  indicate  the  presence 
of  a  definite  shear  stress  concentration  at  the  neutral  axis.  For  all 
beams  the  values  listed  in  TABLE  VI-6A  indicate  that  at  the  neutral  axis 
there  are  no  reasonable  gradients  in  shear  stresses,  between  intact 
sections  in  the  vicinity  of  the  delaminat ions .  Therefore,  the  extent 
of  the  horizontal  shear  stress  concentrations  cannot  be  determined  from 
these  computed  shear  stresses. 

Even  though  the  strain  distribution  in  the  vicinity  of  the  delami¬ 
nation  is  complicated,  it  should  be  possible  to  obtain  valid  shear  stress 
values  from  this  type  of  analysis  on  longitudinal  strains.  In  this 
investigation  a  second  degree  parabola  is  fitted  to  the  recorded  strains, 

/r. 

but  in  very  few  cases  do  these  five  strains  define  a  perfect  parabola. 

A  slight  variation  in  any  one  of  these  five  strains  influences  the 
equation  of  the  best  fit  parabola  to  a  considerable  degree.  Due  to  the 
nature  of  the  computations,  a  small  variation  in  strain  distribution 
produces  a  large  accumulative  error  in  the  shear  stress  at  the  neutral 
axis.  Perhaps  better  results  would  be  obtained  if  deeper  beams  were 
instrumented  and  tested.  With  a  large  number  of  strain  readings  over 
the  depth  of  the  section  a  more  accurate  strain  distribution  might  be 
obtained.  Theoretical  and  experimental  shear  stress  values,  which 
different  computational  methods  indicate  to  exist  in  beam  XY-425,  are 
illustrated  in  FIGURE  6.26.  The  theoretical  shear  stresses  computed  by 
Guyon's  and  Lee's  influence  coefficients  (APPENDIX  C)  indicate  sizable 
shear  stress  concentrations  exist  for  a  distance  of  d/3  (2  inches)  away 
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from  the  end  of  the  de lamination.-  For  greater  distances  these  two 
methods  yield  shear  stresses  which  are  lower  than  the  average  horizontal 
shear  stress,  computed  on  the  undelaminated  neutral  surface.  This  trend 
is  common  for  all  D/L  ratios.  At  a  distance  d/3  from  the  end  of  the 
delamination,  the  two  compatible  values  of  the  shear  stress  at  the 
neutral  axis  are  about  twenty-five  per  cent  greater  than  the  average 
shear  stress.  However,  since  the  occurrence  of  this  compatibility  of 
shear  stresses  at  the  neutral  axis  is  an  exception  rather  than  the  rule, 
no  conclusions  can  be  drawn  with  respect  to  the  magnitude  of  the  stress 
concentration.  In  examining  the  strain  distributions  obtained  for  the 
delaminated  beams  (FIGURES  6.4  to  6.13),  the  extent  of  section  warping 
indicates  shear  stress  concentrations  do  exist  at  sections  situated 
further  than  d/3  away  from  the  end  of  the  delamination.  For  low  D/L 
ratios  warping  extends  to  d/2  (3  inches)  from  the  end  of  the  delamination 
while  for  high  D/L  ratios  warping  extends  to  3/2  d  (9  inches) . 

In  the  delaminated  beams  a  vast  majority  of  sections  exhibit 
outer  tensile  strains  about  ten  per  cent  greater  than  corresponding 
compressive  strains.  This  tendency,  present  at  all  loading  stages, 
suggests  a  slightly  higher  modulus  of  elasticity  in  compression  than 
in  tension. 

At  strains  much  below  compression  yielding,  the  neutral  axis 
in  a  number  of  beams  shifted  down  towards  the  tension  face  of  the  beam. 

As  for  the  control  beam  X-003,  a  reasonable  explanation  of  this  trend 
would  be  a  transformed  section  based  on  an  increased  modulus  of  elasticity 
with  an  increase  in  specific  gravity.  Examination  of  TABLE  VI-5  indicates 
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that  the  majority  of  beams,  in  which  this  shift  occurred  have  denser 
laminations  in  the  lower  portion  of  the  beam  as  compared  to  the  top. 
Therefore,  the  neutral  axis  of  the  transformed  section  would  be 
located  below  the  centroidal  axis  of  the  beam. 

FIGURE  6.20  gives  a  comparison  of  strains  determined  by  the 
2-inch  Demec  gauge  and  SR-4  Type  A-3  strain  gauges.  The  minute  differences 
in  values  obtained  can  be  attributed  to  observation  errors. 

ii .  Strains  in  Delaminated  Sections 

The  method  developed  by  Huggins,  Aplin  and  Palmar  for  determining 

stress  distributions  in  delaminated  sections  was  employed  to  determine 

theoretical  strain  distributions  shown  in  FIGURES  6.14  to  6.18.  A  modulus 

£ 

of  elasticity  of  1.72  x  10  psi,  obtained  from  the  control  beam  tests, 
was  used  in  the  calculations.  The  flexural  stress  distributions,  which 
are  different  for  beams  with  different  D/L  ratios,  must  be  compatible 
with  the  equilibrium  of  the  particular  section  under  consideration. 
Therefore,  the  value  of  the  elastic  modulus  applied  to  the  theoretical 
stress  distribution  will  only  affect  the  magnitude  of  the  theoretical 
strains  and  not  the  actual  distribution  of  these  strains  across  the 
section . 

Beams  XZ-41  and  XY-41,  containing  l'-0"  long  delaminations , 
show  very  good  agreement  between  theoretical  and  experimental  strain 
distributions  (FIGURE  6.14).  For  beams  XZ-42  and  XY-42,  with  2,-0" 
long  delaminat ions ,  the  theoretical  and  experimental  strain  distributions 
are  in  good  agreement  for  sections  1  and  2  (FIGURE  6.15).  At  section  3 
in  both  beams  the  experimental  strain  distribution  is  not  geometrically 
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similar  to  the  theoretical  distribution.  Experimental  strain  distri¬ 
butions  for  beams  XZ-425  and  XY-425  differ  considerably  from  the  theo¬ 
retical  distributions  (FIGURE  6.16).  However,  the  agreement  is  generally 
very  good  for  beams  XZ-43  and  XY-43,  containing  3'-0"  long  delaminations 
(FIGURE  6.17).  FIGURE  6.24A  indicates  the  strains  at  delaminated 
section  3  increase  almost  uniformly  with  increase  in  load.  The  strain 
distributions  for  beams  XZ-44  and  XY-44,  containing  4'-0"  long  delaminat ions , 
agree  with  theoretical  values  at  the  outer  fibers.  Progressing  towards 
the  delamination  plane  this  agreement  begins  to  diverge. 

At  most  delaminated  sections  the  measured  strains  were  of 
lower  magnitude  than  the  theoretical  strains.  This  difference  would  be 
rectified  if  a  higher  modulus  of  elasticity,  as  indicated  by  compression 
and  tension  test  specimens,  is  used. 

An  examination  of  beams  having  the  same  delaminat ion  length, 
but  a  different  type  of  delamination,  indicates  there  is  no  significant 
difference  in  the  strain  distributions  at  delaminated  sections.  This 
implies  that  friction,  in  the  wood  contact  type  of  delamination,  has  no 
apparent  effect  on  the  forces  in  the  delaminated  sections.  In  beams 
having  low  D/L  ratios  agreement  between  experimental  and  theoretical 
strain  distributions  is  generally  good.  However,  for  higher  D/L  ratios 
this  correlation  diverges.  The  initial  assumption  of  plane  sections  re¬ 
maining  plane  after  bending  is  definitely  not  valid  for  sections  at  the 
termination  of  the  delamination.  The  pronounced  warping  of  the  sections 
for  high  D/L  ratios  is  most  likely  responsible  for  the  deviation  in  the 


strain  distributions. 
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For  most  beams  the  correlation  between  experimental  and  theo¬ 
retical  strain  distributions  improved  at  sections  removed  from  the  end 
of  the  delamination.  The  deviation  of  the  theoretical  strains  from 
the  actual  strains  near  the  end  of  the  delamination  affects  the  magni¬ 
tude  of  the  shear  stress  concentrations,  computed  using  Guyon's  or  Lee's 
influence  lines. 
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CHAPTER  VIII 


CONCLUSIONS  AND  RECOMMENDATIONS 


8-1  Conclusions 


The  most  significant  conclusions  based  on  the  results  of  this 
investigation  are  summarized  as  follows: 

1.  Non-linear  longitudinal  strain  distributions  at  intact  sections 
of  delaminated  beams  indicate  shear  stress  concentrations 
exist  in  the  vicinity  of  a  delamination.  The  extent  of  this 
non-linearity  in  strain  is  more  pronounced  in  beams  having 
large  delamination  length  to  beam  length  ratios. 

2.  It  is  not  possible  to  accurately  compute  shear  stresses 
from  equilibrium  considerations  using  longitudinal  strain 
distributions  obtained  in  this  investigation.  The  relatively 
small  number  of  strains  recorded  across  the  beam  section  do 
not  satisfactorily  define  the  strain  distribution. 

3.  The  theoretical  and  experimental  strain  distributions  at 
delaminated  sections  agree  reasonably  well  for  low  delami¬ 
nation  to  beam  length  ratios.  However,  this  correlation 
decreases  with  increasing  delamination  length  to  beam  length 
ratios . 

4.  At  delaminated  sections  it  appears  that  there  is  no  appreciable 
friction  between  the  contacting  delamination  surfaces. 
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8-2  Recommendations 

The  results  of  this  investigation  indicate  that  tests  on  deeper 
beams  may  be  warranted.  An  increased  section  depth  may  facilitate  a  more 
accurate  determination  of  the  longitudinal  strain  distribution  across 
the  section. 

The  possibility  of  developing  further  theoretical  shear  stress 
distributions  from  non-linear  longitudinal  strain  distributions  should  be 
investigated . 
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APPENDIX  A 


DETERMINATION  OF  MOISTURE  CONTENT  AND  SPECIFIC  GRAVITY 

Specimens  were  weighed  prior  to  being  dried  in  an  oven  at  100°C 
for  36  hours.  Upon  removal  from  the  oven  the  specimens  were  placed  in 
a  desiccation  cabinet  for  2  hours  and  then  weighed.  The  volume  of  the 
specimen  was  obtained  by  a  mercury  displacement  method.  A  container 
was  filled  with  mercury  until  the  volume  of  mercury  equaled  the  exact 
volume  of  the  container.  The  specimen  was  then  immersed  in  the  mercury 
and  the  mercury  which  overflowed  from  the  container  was  weighed.  Re¬ 
peating  this  procedure  for  any  one  specimen  indicated  the  maximum  vo¬ 
lumetric  deviation  of  displaced  mercury  specimen  was  one  percent.  Moist 
weight,  dry  weight  and  weight  of  displaced  mercury  were  read  into  the 
IBM  7040  Digital  Computer  and  the  moisture  content  and  specific  gravity 
were  calculated  by  the  program  on  page  A2. 
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$  I B SYS 

SJOB  180027  E.  F.  COLLINS  PROGRAM  ONE 

ST  IMF  2*100  0 

$  I B JOB  REDUCE 
$ IBFTC  RESULTS  NODECK 

C  COMPUTATION  OF  moisture  content  and  SPECIFIC  GRAVITY 
COMMON  MW( 8* 1 5 ) *DW( «  *  1 5 ) ♦ WM ( 8 ,l 5 ) *SG (8*15 ) *MC( 8  *1 5 ) 
INTEGER  B 
REAL  MW*MC 
DO  500  8=1*13 

RE AD ( 5  *  1 00 )  ( MW ( L ♦ B ) *DW(L*B)  *WM(L*B)  *  L  =  1  *  4 ) 

RE AD ( 5  *  1 00 )  (MW(L*B) *DW(L*B) *WM(L*B)  *L  =  5*8) 

WR  I  TF ( 6  *  20 1 ) 

DO  501  L= 1  *  8 

779  MC(L*B)=(MW(L*B)-DW(L,B) )/DW(L*B) 

778  SG(L*B)=DW(L*B)*13.55/WM(L,p) 

501  WRITE(6*200)L*MW(L,B) *DW(L,B) *WM(L,B) *MC(L,B) *SG(L,B) 

500  CONTINUE 

201  FORMAT ( 5HLBEAM ) 

200  FORMAT! IX* I5*5F10.3) 

1 00  FORMAT ( 1 X, 2F4. 1  *F6. 1 *2FA.l  *F6. 1  *2F4, 1  * F6 • 1  *  2 F4 . 1  * F6 . 1  ) 
STOP 
END 

SENTRY  RESULTS 


. 

eve  8u 

i  A&doaq  .3  * j 

32U03fl  80L3 1  * 
5T3C10&1'  2TJU339  DT'I8!$ 

8  *3t>3TM1 

i  ;"  •  \ • 

l'  T  *  f*S  0  0  r*  (  ) 

*  '  > »  <  •-»'••••  f  ■  /(»(  •■•  )  ■  ■;  ft'’}  o  3  ■ 

( ics.anM  fnv 

•  I  J  n 

f  •  ♦  J)V3\  ,  JSlvMl-  (P;,  J)  -)M  ?rr 

<  Q  *  J )  Mw  \  5  ?  .  C  1*  <  8  *  J  >  WO  -  {  g ,  J  )  ap,  gy  y 
n  J  r  .  <g>  !>  '  ’  roo:<<3  ~it  P\  f 

Bui-irn  ooe 

'  i  N1  A  3  8  J  H  )  fOp 

<  _  •  r* r  ‘  )  •'  • 

•  f  •  ’>  1  f  I  *A3S  #  f  ♦  <&  t«  f  #.\ls  t  r  r  *  r»£l*  f  *Xf)  TAWflO  h  00 ; 

qoTe 

erJ»>23#  YP^a* 


APPENDIX  B 


COMPUTER  PROGRAM  FOR  SHEAR  STRESSES 

The  reduction  of  strain  gauge  readings  and  the  subsequent 
computations  in  determining  shear  stresses  are  time  consuming  and  re- 
petative.  Therefore,  solutions  were  obtained  by  means  of  a  computer 
program.  This  program  was  written  in  Fortran  IV  for  the  IBM  7040  Digital 
Computer.  A  very  simplified  flow-chart  and  the  actual  program  as  shown 
on  pages  B 2  to  B8 . 
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FLOW  CHART  FOR  SHEAR  STRESS  COMPUTER  PROGRAM 
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S  IBSYS 

SJOB  180027  E.  F.  COLLINS  PROGRAM  TWO 

STIME  8*50QO 

S I B JOB  REDUCE 
SIBFTC  RESULTS  NODECK 

C  COMPUTATION  OF  SHEARING  STRESSES  FROM  EQUILIBRIUM  CONSIDERATIONS 

COMMON  P( 1 2 )  *DR ( 10* 12  *2 ) *ST ( 1 0 , 12  *2 )  *NA ( 1 2  *2  >  *NA2  (12*2)* 

1  X  (  10*12*2) *  Y ( 10*12*2)  *YS(63*12*2) *YST ( 63  *  1 2 ♦ 2) ♦ XB ( 63 ♦ 1 2  ♦  2  )  . 
2F(63*12*2)*M (63*12*2)  *XBC(63*12*2)*Y SC (63*12*2)  * VSTC ( 63  *  1 2 ♦ 2 )  * 
3FC(63* 12 *2)* MC (63*12*2)  *DFC (63*12*2)  * YVC ( 63  *  1 2  *  2 ) *SVC (12*2). 
4NES.NEF  *NP*N1 1  * Nl 2 1  * N2 22  * N24 1  * E *  EC  * ET * D * B * NX  * L , XS ♦ XSL * LO , Nl 0  * DNA . 
5C*GL  *TL*GLX,H*NGA1 *NGA  *  XX  *  XXX ♦ AMR ,BMR , CMR , DC  * NNA ( 12 ) 

REAL  P*NX,NA*NA2*NF*NGA*NF2*NF3 
INTEGER  XS*GL,XSL*H,X2*C 
READ ( 5  *  1 0 1 )  N0 
READ(5*105)EC*ET 
READ ( 5  *  1 06 ) D  *B 
READ(5*101 )NP 
DO  200  NBC  =  1  *  NB 
WR I TE ( 6  *  1 09  )  NBC 
109  FORMAT ( IX*  17) 

RE  AD ( 5  *  1 00 )  NES.NEF 
NX  =  3  0 • 

DC  =  0*10 
TL  =D/8 • 0 
N 1 1  =  1 
N24 1 =62 
L=  1 

DO  201  XS  =  NES  *  NEF 

201  READ (5. 103) (DR(GL,XS.L) *GL=1*9) 

DO  202  LC  =  2  *  NP 
WRITE(6*102) 

L  =  2 

DO  300  XS=NES*NEF 

300  RE AD (5*103) (DR(GL*XS*L) *GL=1 *9) 

DO  203  XS  =  NES  *  NEF 
LO=  1 

DO  204  GL=1 *9 

204  ST ( GL  *  XS*  L )  =  ( DR ( GL  *  XS  *L )-DR ( GL  *XS  *LO ) ) *0. 03  8 
N4  =  4 
N5  =  5 
N6  =  6 
Nl 0=1 0 

S  T  4= ( ABS( ST ( N4  *  XS  *  L )  )  ) 

ST  5= ( ABS ( ST(N5*XS*L)  )  ) 

ST6=(ABS(ST (N6*XS*L)  )  ) 

I F ( ST ( NS ♦ XS * L ) . LE . 0* 00 )  GO  TO  253 
DNA=-TL*ST5/ (5T5+ST4 ) 

GO  TO  281 

253  DNA=TL*ST5/(ST5+.ST6) 

281  DO  254  GL  =  1 *9 

254  ST ( GL  *  XS*  L ) =ABS (ST ( 6L  *  XS  *  L )  ) 

NA ( XS  *  L ) =D/2*0+DNA 
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■  N A  2 ( XS  »L )  =  D-NA ( XS*L  ) 

GL  X  =  NA ( X5  *  L ) +TL 

NF=NA!XS.L)/TL+1.0 
NI  =NF 

DO  282  GLsl.Nl 
6LX=GLX-TL 
X ( GL  ♦  XS  *L ) =GLX 
PPP  Y(GL*XS*L)=ST ( GL  *  XS  *  L  ) 

H=  1 

NGA I =N I +1 
N 1  0=  N  I  + 1 
X! N10,XS*L)=0.0 
Y( N10.XS.L )=0.0 
NGA=NGA I 
CALL  SRP1 
CALL  SRP2 
CALL  SRP3 
NF2  =  9 • O-NF 
N I 2=N I +1 
N I 3=NF2 
NF3  =  N 1 3 

GLX=NA2(XS»L)-NF3*TL 

DD2=D*0,5 

DD3*NA ( XS  *  L ) 

DD4=ABS(DD2-D03 ) 

I F ( DD4*  EQ*  0*  00 )  GLX  =  0.73 
NG A  I  =  1 1 -N I  2 
NG A=NGA I 
N9  =  l 

DO  283  GL  =  Nl 2*9 
N9=N9+1 

X( N9,XS»L)=GLX 
Y(N9*XS*L)=ST(GL*XS*L) 

283  GL  X=GLX+TL 
H=  2 
N  1  0  =  1 

X ( N10*XS*L )=0*0 
Y(N10*XS*L)=0#0 
CALL  SR PI 
CALL  SRP2 
CALL  SRP3 
203  CONTINUE 
CALL  SRP4 
DO  303  C=N 1 1  *  N241 

303  WRITE (6* 104) (XBC(C*XS,L) *YVC(C*XS*L) ,XS=2*NEF) 

202  CONTINUE 
200  CONTINUE 

100  FORMAT! IX *2 13) 

101  FORMAT! IX*  13) 

102  FORMAT ( 1HK  ) 

103  FORMAT! 1X*9F7.3> 

1 04  FORMAT ( lX*F11.2*F7.1*Fll*2*F7.1*F11.2*E7.1*F11.2*F7.1»Fll.2f 
1  FI  1.2*F7. 1 ) 
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105  FORMAT! 1X*?F8.0) 

106  FORMAT ( 1 X  *  2F 5 • 2 ) 

1 07  FORMAT ( IX* 10F6.0) 

STOP 

END 

SI6FTC  ABCD1  NODECK 
SUBROUTINE  SRP1 

COMMON  P (  1 2  )  »DR ( 10 » 1 2 » 2 ) • ST ( 1 0 » 12  »2) »NA ( 12 *  2) * NA2  ( 12*  2 ) » 

1 X ( 10*12*2) *Y(10*12*2) »YS(63*12*2) *YST(63*12*2) *XB(63*12*2) * 

2F( 63*12*2 ) *M( 6  3 *12*2 ) *XBC( 63* 12*2) * YSC ( 63  *  1 2 , 2 ) * YSTC( 63 , 1 2  * 2 ) ♦ 

3FC  ( 63 *12*2)* MC (63*12*2)  *DFC (63*12*2)  * YVC ( 6 3  *  1 2  *  2 ) *SVC( 12*2)* 
4NES*NEF*NP,N11 ,N121 *N222*N241  * E *  EC , E T , D * B *  NX , L , XS * XSL , LO * N 1 0 . DNA * 
5C  *  GL  *  TL *GLX , H , NGA I  ,NGA , XX , XXX  *  AMR ,BMR ,CMR ,DC *NNA ( 12) 

REAL  AMR ,BMR,CMR,M1 ,M2  *M3  *M4 ,M5 ♦ M6 * M7 *M8  *DD * EE *BB , CC , X ♦ Y , NGA 

integer  G*C*L*XS 

Ml =  0  •  Q 

M2  =  0 • 0 

M3  =  0 • 0 

M4  =  0 • 0 

M5  =  0 • 0 

M6  =  0  •  0 

M7  =  0  •  0 

M8=0.0 

DO  600  G= 1  *  NGA I 
M1=M1+Y(G*XS*L) 

M2=M2+X(G*XS*L) 

M3=M3+( (ABS( Y(G*XS*L) ) ) **2 • 0 ) 

M4=M4+( (ABS(X(G*XS,L) ) )**2.0) 

M5=M5+( (ABS(X(G*XS*L) ) )**3-0) 

M6  =  M6-M  (ABS(X(G*XS*L)) )**4.0) 

M7=M7+(X(G*XS*L)*Y(G*XS*L) ) 

600  M8  =  M8+(  (ABS(X(G*XS*L)  )  )  **2 • 0 ) * ( Y ( G * XS ♦ L )  ) 

DD= ( NGA* ( M4*M£-M5*M5 ) ) - ( M2* ( M2*M6~M5*M4 )  )  +  ( M4* ( M2*M5~M4*M4 ) ) 

EE  = (  Ml  *( M4*M6-M5*M5 )  )  - ( M7* (M2*M6”M5*M4 )  )  +  ( Ms* ( M2*M5“M4*M4 )  ) 
BB=(NGA*(M7*M6~M8*M5) )- (M2* ( M] *M£-Mq*M4 ) ) + ( M4* ( Mi*M5~M7*M4 ) ) 

CC= ( NGA*( M4*M8"M5*M7 ) ) - ( M2* (M2*M8"M5*Ml ) ) + ( M4* ( M2*M.7"M4*Ml ) ) 

AMR=EE/DD 

BMR=BB/DD 

CMR=CC/DD 

RETURN 

END 

SIBFTC  ABCD2  NODECK 
SUBROUTINE  SRP2 

COMMON  P ( 1 2 )  *  DR (10*12*2) *ST(10*12*2)  *NA(12*2)*NA2(12*2)  ♦ 

1 X (  10*12*2)  *Y( 10*12*2)  *YS(63*12*2) *YST (67*12*2) *  XB (67*12*2)  ♦ 

2  F ( 63*12*2)  *M ( 67  *  1 2  *2  )  *X PC (63*12*2) *YSC(63*12*2)  *YSTC(63*12*2) ♦ 
3FC(63*12*2) *MC(63*12*2) *DFC(67*12*2) *YVC(63*12*2) ♦ S VC ( 1 2 ♦ 2 ) ♦ 

4NES  *NEF  *NP  >N1 1 *Nl 2 1  *  N222*N24 1  *  E*  EC*ET  *  D  *B  *NX ♦ L , XS  *  XSL  *LO  *Nl 0  *  DMA  » 
5C*GL ,TL*GLX,H*NGAI  ♦  NGA *  XX , X X X , AMR  * BMR , CMR , DC  * NN A ( 1 2 ) 

REAL  NX*M*MI *ME*NA*NA2 
INTEGER  H, XS*XX*XXX* XI *C*XX1 
XI  =63 

IF  ( H#  EQ«  2 ) GO  TO  602 
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Cl =NA(XS,L)/DC 
NL  =  C1 
C3  =  N 1 
N 1 Z=N 1 1+1 

XX  =  N  1 1 
XXX=Nl+2 
N1  21 =  XXX 
NNA ( XS) =N1 21 

XXI  =  0 
E=  EC 

GO  TO  603 

602  XX1=XXX+1 
N22  2  =  XX 1 
XX=N222 
XXX=N241 
C2=N1+1 

XB ( X 1  *  XS  *  L ) =C2*DC-NA ( XS . L ) -DC 
E=ET 

603  YB( X1,XS*L)=0.0 
YST(X1,,XS,L)=0.0 
M ( XI ,XS*L)=0.0 
00  601  c=xx*xxx 

IF (C.EQ.N12)  GO  TO  608 
IF (C.FQ.XX1 )  GO  TO  611 
IF (C.E0.N1 1 )  GO  TO  636 
XB( C ,XS  *L ) =XB( XI ,XS*L )+DC 
DC  1=1.0 
GO  TO  610 

608  XB(C*XS*L)=NA(XS*L)-C3*DC 
DC1=ABS(XB (C*XS*L) ) /DC 
GO  TO  610 

611  XB(C,XS*L)=C2*DC-NA(XS*L) 

DC1=ABS(XB(C*XS*L)  )/DC 

610  YS(OXS>L)=cMR*( ( ABS( XB(C,XS*L) ) )**2.0)+BMR*xB(C,XS»L)+AMR 
YST(C,XS*L)=YS(C*XS*L )*E 

F(C,XS,L)=  (YST(C,XS»U+YST(X1  *XS*L)  )  *0.5*DC*B*DCl 
M(C*XS*L)=F(C*XS*L)*( XB ( C , XS * L ) -0 .5*DC*DC 1 )+M(Xl ,XS,L) 

GO  TO  601 
636  F(C*XS*L)=0.0 
M(C*XS*L)=0.0 
YS ( C*XS*L ) =0.0 
YST(C,XS*L)=0.0 
XB(C  ,XS*L) =0.0 
601  X1=C 
RETURN 
END 

STgFTC  A0CD3  NODECK 
SUBROUTINE  SRP3 

COMMON  P  (  1  2 )  ♦DR(10»12*2)*ST(10*12*2)  *  NA ( 1 2  *  2 ) *NA2( 12*  2)  ♦ 

1X(  10* 1 2 *2) *Y( 10*12*2) *YS( 63 *12*2) *YST( 63 *12*2) *XB (63*12*2). 

2F( 63*12*2) *M( 63 *12*2) *X BC (63*12*2) *Y SC (63 *12*2) »YSTC( 63 *12*2)* 

3FC (63*12*2) * MC (63*12*2)  * DFC ( 6 3 ♦ 1 2  *  2 )  ♦ YVC ( 63  *  1 2  *  2 ) ♦ SVC ( 12*2) * 
4NES*NEF ,NP  »N1 1 *N121 > N222 * N2 4 1  *  E*  EC, FT  * D ♦ B *  NX , L , XS * XSL  * LO , N 10 * DNA , 
5C,GL,TL*6LX,H,N6AI * MGA *XX , XXX , AMR , 0MR , CMR ,DC , NNA ( 12 ) 

REAL  NA 
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INTEGER  XX»XXX*X2*K*X5*L*H 
I F ( H.EQ.2 )  GO  TO  850 
X  X  =N  1  1 
21=1 .0 
XX  X  =  N1 21 
X2  =  N1 21 

<  =  -l 

GO  TO  852 

850  XX=N222 
XX  X  =  N24 1 
Z= -1 .0 

K  =  1 
X2=XX 
GO  TO  832 

852  DO  851  NN=XX*XXX 

X0C(NN*XS*L)=NA(XS*L )-XB(  X2*XS*t  )*L 
VSCt NN,XS,L)=Y$( X2*XS*L) 

YSTC(NN,X$*L)=YST(X2*XS*U 
FC (NN,XS*L)=F(X2*XS*L) 

MC (NN,XS*L)=M( X2  *  XS  ♦  L. ) 

851  X2  =  X2  +  K 
RETURN 
ENO 

SIBFTC  ABCD4  NODECK 
SUBROUTINE  SRP4 

COMMON  P(1 2) * DR ( 10*12*2) *ST ( 10,12*2) *NA( 12*2)  *NA2(1  2*2)  * 

1X( 10*12*2) *Y(1 0*12*2)* YS( 63 *12*2) * YST ( 63  *  1 2 , 2 ) * XB ( 63* 1 2 *2 ) * 

2F  (  63*1  2*2.)  *M(  63*1  2*2)  *X»C(  63*1  2*2)  *  YSC  (  63  *  1  2  *2  )  *YSTC(  63*  1  2  *2  )  * 
3FC  (63*1  2*2)  *MC(  63*  12*2)  *DFC  (  63  *  1  2  *  2  )  *  YVC  (  63  » 1 2  *  2  )  *SVC  ( 12*2  )  * 
4NES*NEF ,NP ,N11 *N121 *N222*N241 *E*EC*ET*D,B,Nx*L,XS,XSL,LO,N10,DNA 
5C  *GL ,TL*GLX,H*NGAI *NGA *  XX  *  XXX , AMR ,BMR * CMR ,DC *NNA  (  1 2 ) 

INTEGER  XS1,XS*C*CC 
DO  704  XS=2*NEF 

xsi=xs-i 

K 1 =N 1 1 
K2=NNA( XS ) 

K  3  =  K 1 

<  4  =  <  2 
<8  =  1 

K  6  =  - 1 
GO  TO  719 
719  C=Kl+<6 

DO  702  CC=K3*K4 
C  =  C  +  K5 

IF ( C.E0.K1 )  GO  TO  709 
IF (C.EO.K2)  GO  TO  71  0 
GO  TO  713 

709  DFC(C,XS*L)=(FC(C,XS1 , L ) “FC ( C , XS ♦ U ) / ( 7. 0*B ) 

YVC ( C  *  X S ♦ L ) =0.0 

SVC(  XS*L)  =  DFC(C*XS*L ) 

GO  TO  702 

710  YVC(C,XS*L)=SVC(XS*L ) 

GO  TO  702 

713  DFC(C*XS*L)=(FC<C,X$1 ,L)-FC(C,XS*L) )/(2#0*B) 

YVC( C*XS»L ) =5VC( XS*L ) 
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SVC( XS  » L )  =  DFC( C*  XS  *  L ) +SVC ( XS  *  L ) 
702  CONTINUE 

IF(K1.FQ.N11 )  GO  TO  720 
GO  TO  704 
720  K 1  =  N24 1 
K2  =0 

K3=NNA( XS) 

K4=N241 
K  5  =  - 1 
K6  =  1 

GO  TO  71 Q 

704  continue: 

RETURN 

FNO 
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APPENDIX  C 


CALCULATION  OF  THEORETICAL  SHEAR  STRESS  CONCENTRATIONS 


The  calculations  of  the  theoretical  shear  stress  concentrations 
(FIGURE  6.26),  based  on  Guyon's  and  Lee's  methods,  are  outlined  in 
this  appendix.  The  following  sample  calculations  illustrate  the 
determination  of  the  horizontal  shear  stress,  based  on  Lee's  method, 
at  a  distance  d/6  (1/2  inch)  from  the  end  of  the  delamination  in 
beam  XY-425. 

The  distribution  of  flexural  stresses  across  a  section  at 
the  end  of  the  delamination  must  be  obtained  by  the  simultaneous 
solution  of  Equations  1.1  and  1.4  as  follows: 

Rt  (a  +  z)  -  %  -  Mb  -  y  =  0  d-D 


x  =  2 

x  =  0 


(1.4) 


where  a  =  21.0  in. 

a^  =  1.0  in. 

z  =  15.0  in. 
b  =  3.0  in. 
d  »  6.0  in. 
k  =  0.5 


AT  =  AB  =  9.0  in. 


A  =  18.0  in. 


I  =  54.0  in. 


IT  =  Ig  =  6.75  in 


=  2500  lbs. 
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Since 


k  =  0.5 

M  =  M 
B  T 

C  =  T 


=  0.5  R 

L 


Substituting  these  values  into  Equations  1.1  and  1.4: 

110  -  2Mj.  -  3C  =  0  (1.1A) 

20.8  -  2 . 22Hj.  +  1. 11C  =  0  (1.4A) 


Solving: 


C  -  T  =  17.8  kips 

Hp  =  Mg  =  18.25  in-kips 


At  the  end  of  the  delamination  where  x  =  15  in.: 


Mtx  =  Mj.  -  VT  X 

MpX  =  MgX  =  -0.50  in-kips 

The  flexural  stresses  in  the  top  and  bottom  delaminated  portions  are 


f  =  + 
max .  — 


41-j- 

1,987  psi  at  the  delamination 


fmin  =  +  1,967  psi  at  the  outer  fibers 
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The  average  horizontal  shear  stress  on  the  undelaminated  surface  is: 

C 

vave  =  b(a]L  +  a)  =  270  Psl 


In  order  to  compute  the  shear  stress  concentration  on  the 
neutral  surface  at  d/6  from  the  end  of  the  delamination,  each  delami 
nated  beam  portion  is  divided  into  ten  elements.  Influence 
coefficients  (FIGURE  1.7)  are  obtained  for  each  element  knowing  its 
distance  from  the  neutral  axis.  The  shear  stress  at  the  neutral 
axis  due  to  the  influence  of  a  single  element  is: 


where 


0.05  bd  fnCn 


vn  = 


bd 


0.05  fnC 

’n 

n  " 

element 

"n 

element 

"n 

'n  = 

element 

"n 

shear  stress 
fiber  stress 
influence  coefficient 


For  the  twenty  elements  through  the  depth  of  the  beam  the 
total  shear  stress  is: 

n  =  20 

v  =  0.05  £  fnCn 

n  =  1 

Since  the  computations  for  the  two  delaminated  portions  are  identical 

n  =  10 

v  =  0 . 10  £  fnCn 

n  =  1 

The  computations  are  tabulated  below  in  TABLE  Cl. 


: 


t:x  eai  ?:  j  aril  £»ai/q  >-  o3  yabio  nl 

■_r:.s  ,00 Jtj2.fi,;  :.(  >.>  ul  h  ln  rH  ro  ri  d\b  3ft  9»-.'  tue  X  .?3  »r 

, : .  '  ivoirt  J  ■.;  :.v  a  rioas  *  :  ■  ?. - u  iBlxfo  «ia  ("  i  *.*i:  rri)  ,rr-.  :  iilon 

Hifsvfj  9:i3  ift  ae^i3t  i8»j-  h  ®;?T  .a£Xft  .  »*  ^on  ar  1  mot'i.  wilalb 
: #1  jnarafd  ,  nie  i  }o  ar>nat>.  ni  vrta  o3  *ub  ntxB 

,7on5  bcf  ro  0 


bi 


*  nv 


'  n  n 


0  ,1  eo.o  - 


•  J'Je  ri  *.  M.n  '  '  •/!  il  9  ~  n>' 

...  .  a  -■'*d  r,  1  u  iJ  3  -  ,  5 


ifia  i  o  /  1  it 900  o  >n^u  i  ill  V’  inoan  I  &  ~ 


9?*)rfw 


r,  -  3  '  "0  ')  91  -  SlOSUf*rp  7lf!9W.1  9/tJ  JO" 

:ei  ziiQsjz  n>*rf8  XbJoJ 

f 

0£  -  a 


r. 


3  r.o.o  -  v 

1*3 


X6i  .!•  ir  s  o no i 3  icq.  b^lfti  rirmlsb  ow'  ♦rt.;  10 i  *  k  1 4j:  .luc-moo  orto  ooriJ-’. 


OX  ~  n 

s  01 .0  -  V 
1  =  n 


•  I"'  XJfffAT  nl  wolsd  badr  JucfB3  91B  &rx>X JB Juqwoo  oriT 


TABLE  Cl 


C4 


LEE'S  SHEAR  STRESS  CONCENTRATION  AT  d/6 


Element 

f  (psi) 
n 

C 

n 

f  C 
n  n 

1 

1968 

i 

o 

• 

u> 

o 

-  590 

2 

1970 

-0.35 

-  689 

3 

1972 

-0.20 

-  394 

4 

1974 

i 

o 

H-* 

o 

-  197 

5 

1976 

0.05 

99 

6 

1978 

0.18 

356 

7 

1980 

0.45 

890 

8 

1982 

0.93 

1843 

9 

1984 

2.05 

4067 

10 

1986 

1.02 

2035 

n  =  10 

Z  f  C  =  7420 

n  n 

n  =  1 

v,  „  =  0.10  x  7420  =  742  psi 
d/6 

Similar  computations  were  performed  for  other  sections  and  the  results 
recorded  in  FIGURE  6.25. 
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